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Sumroary 

A theoretical  investigation  is  carried  out  on  the  buckling  of 
sandwich  cylinders  under  combined  compression,  torsion,  and  bending 
loads.  The  governing  differential  equation  is  solved  by  using 
Galerkin’s  method.  The  interrelationship  obtained  between  the  criti- 
cal loads  is  plotted  in  the  form  of  non-dimensional  interaction  curves. 
In  the  limiting  cases  of  axial  compression  alone,  torsion  alone,  bend- 
ing alone,  and  combined  bending  and  axial  compression,  the  results  agree 

1-3 

with  those  obtained  previously. 
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Symbols  end  Unite 

radius  of  cylinder  to  mid-plane,  in0 
hG,  lbs0  per  in. 

flexural  stiffness  of  isotropic  sandwich  cylinder,  in.  lbs. 

D « Eth2/2(1  - v2) 

Young's  rcoaulus  for  the  face  material  lbs.  per  sq.  in. 

2tE/a2 

shear  modulus  of  the  core  material, lbs.  per  sq.  in. 
depth  of  isotropic  sandwich  plate  measured  between  middle  planes 
of  faces,  in. 

moment  of  inertia  of  the  sandwich  cylinder  about  its  diameter 

length  of  cylinder,  ln0 

bending  moment,  in.  lbs. 

number  of  half  waves  longitudinally 

resultant  normal  forces  in  x-  and  y-directions,  lbs.  per  in. 
resultant  shearing  force,  lbs.  per  in. 


applied  normal  force  in  axial  direction,  lbs.  per  in. 


number  of  full  waves  circumferentially 
torsional  moment,  in0  lbs0 
thickness  of  the  faces,  in. 

displacements  in  x-,  y-,  z-direction3,  respectively,  of  a point  in 
middle  surface  of  cylinder,  in. 
rectangular  coordinates  (Fig.  1) 

Poisson's  ratio  fur  the  face  material 


circumferential  coordinate 

a2  a2  2 a4 

operator.  ( — -x  + — x ) - (=-r  +2  - ? 

3x  5y  3x4  8x  dy 


*4  ^ a2  a2  2 

+ *■ 

8y4  8x  a^a© 


.a2  a2  v4 

operator,  v — j + — “ ) 

ax  ay* 


/8 

W 
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Introduction 
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In  previous  papers  the  senior  author  and  his  associates  have 

investigated  the  buckling  behavior  of  sandwich  cylinders  under  axial 
compression,  torsion,  bending  and  combined  bending  and  axial  compression. 
The  investigation  is  now  extended  to  the  buckling  of  sandwich  cylinders 
inder  combined  compression,  torsion  and  bending  loads.  By  solving 
Donnell's  equation  modified  to  include  the  effects  of  transverse  shear 
for  sandwich  curved  plates  and  shells  using  Galerkin's  method,  the  buck- 
ling loads  are  calculated  and  are  plotted  in  the  form  of  interaction 
curves.  In  the  limiting  cases  of  bending  alone,  axial  compression  alone, 
torsion  alone,  and  combined  bending  and  axial  compression,  the  results 
agree  with  those  obtained  previously. 


Formulation  of  The  Problem 
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By  assuming  isotropic  core  material  and  neglecting  the  bending 
rigidity  of  the  faces  about  their  own  middle  surfaces*  the  equilibri- 
um equations  for  element  of  a sandwich  cylinder  have  been  derived  in 
reference  4.  By  neglecting  terms  that  were  regarded  &j  small  by  Donnell*5, 
these  equations  can  be  reduced  to  a single  equation  in  terms  of  the  lateral 
deflection  w only  • - 


Dv8v  + (1  - - V2) 
C 


2tE  84v 
a2  8x4 


(N  “ + N + 2N  — ) 

X 8x2  y 8y2  ^ 8x8y 


(1) 


This  is  the  so-called  Donnell’s  equation  modified  to  include  the  effects 

of  transverse  shear  for  sandwich  curved  plate3  and  shells  first  obtained 

7 

by  Stein  and  Mayers  , 

In  the  case  of  sandwich  cylinders  under  combined  compression,  torsion 
and  bending  loads,  we  have 


Nx  - - + (2Mat/l)  cos  0, 


(2) 


where  is  the  force  per  unit  length  due  to  axial  compression,  M is  the 

bending  moment,  I is  the  moment  of  inertia  of  the  cylinder,  and  T is  the 
torsional  moment.  With  such  a loading,  equation  (l)  thus  becomes 


D v8'-'  ♦ (l  - fi  V2) 

r. 


% 

) 2tE 

e4w  4 

- 

( 

• N,  + 2—  cos  0) 

82w  T 

e2w"|i 

k 2 

~4  "V 
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I a 

dx 

I 

dx  fra 

(3) 


where  y is  taken  as  a Q, 


Galerkla*a  Method 

Equation  (3)  may  be  solved  by  means  of  Galerkin’s  method  as  follows* 
we  first  assume  the  deflection  v of  the  cylinder  after  buckling  in  the 
form  of  a series  that  satisfies  the  boundary  conditions  but  with  unde- 
termined parameters.  For  long  cylinders,  however,  the  boundary  condi- 
tions at  the  two  ends  become  unimportant  and  we  may  assume  the  deflection 
series  without  any  regard  for  the  end  conditions.  For  a cylinder  under 
combined  axial  compression  and  bending,  w may  be  assumed  in  the  following 
forra^, 

w « sin  Z A cos  n Q (4) 

1 n-0  n 

and  for  a cylinder  under  torsion  only,  the  deflection  after  buckling 
is  of  the  following  form*, 

w ■ Z Br  sin  - nO) 

n-1 

■ZB  (sin  cos  - cos  sin  nO)  (5) 

n-1  1 

Guided  by  these  expressions,  ve  shall  assume  the  deflection  of  the  cy- 
linder in  the  case  of  combined  loading  as  follows: 

_ nl  __  nl 

w - sin  —t—  Z A cos  nQ  + cos  —r—  Z 3 sin  n£.  (6) 

1 n-0  n 1 n-1 

With  v expressed  in  a proper  series,  we  shall  nexx  substitute  this  series 
into  (3).  If  the  expression  (6)  happens  to  be  the  exact  solution  of 
equation  (3),  after  substitution  equation  (3)  will  be  identically  equal 
to  zero.  In  general,  this  will  not  be  so  and  the  resulting  expression 
will  be  a function  of  x and  Q which  we  shall  denote  by  Q.  Galerkin*s 
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oquations  for  the  determination  of  the  coefficients  and  ore 
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Q sin  cos  nO  ad©  dx  = 0 


Q cos  sin  n 0 ad©  dx  * 0 


(7) 


(8) 


which,  when  written  out,  become 

l 


♦ *, 


Zf'  f (D  V^w  ♦ p - — y^  2-y  ) gin  2EZ  cos  nQ  fa 

J J ex4  C dx4  1 

0 

f ( y4  ) 3in  cos  n©  ad©  dx 

7 Ox*  C 8x*  1 


0 

2ir 


0 0 


2M 

I 


( y4  ) sin  y — cos  n©  cos  © ad©  dx 

Ox  C 6x*  L 


Tra 


4 82v  _ D ^6  8^v 

6x00  C V 6x6© 


J 

I 


S 

) sin  cos  n©  ad©  dx  * 0, 


(9) 


in  which  F *•  2tE/a  , and 


» - 


4 

* 


(S 


sin  nO  ad©  dx 


DF  2 6^w  ->  bttx 
C Ox*  1 


norx 


n©  ad©  dx 


ex* 


y”  *-*  } cos  ~r“  sin  nO  cos  © ad©  dx 
Ox 


2 

( V74  ^ ) cos  sin  n©  ad©  dx  • 0 . 


(1C) 


Substituting  the  assigned  function  v into  these  equations  and 
carrying  out  tbs  integration,  vs  obtain 


Bi  X 2 ; n2)4  , Tjt  . PF  X4  (t?  j n2) 

8 4 _ 6 
a a Ca 

2/  -\2  . _2\2  «»  n \ 


N.  X ( X ♦ n ) 


Ca 


Tn  X(  X2  + n)Z  . Wn\(\‘ 
3 . 10 

irCa 


^2  . _2\2 
,2^j2  . _2X2 

la 


Mi  * p 1 * *n}  (vi  * w - o. 

(U) 


ICa 


A2  ♦ 

I 

CF\ 

No 

1 

J 

1 1 

B 

n 

^2' 

^V-l 

and 


PIJLjjlI4  , ol  , PF  \4i  A2  « a2) 
§ 4 Tb 


Ca 


■»  2/  -.2  2%2  v » \ 2/  2 2%3 
X (X  * n ) N^D  X ( X + n ) 


Ca 
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B_ 


+-b!£  + mlLX 


ira 


uCa' 
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u££ 


Mat 


>2(  A2  » nV  , M»tD  Xa(^  ♦ nV  j (g^ 


* W ‘ °* 


(12) 


Instead  of  working  with  these  two  equations.  It  Is  found  oonren- 
ient  to  combine  them  into  the  following  ones*  Thus,  by  subtracting 
equation  (12)  from  (11),  we  obtain 

D(  X2  ♦ »2)4  ,'Tjt  , nr  \4<  X2  « »2) 

a a*  Ca 

X 2(  X2  ♦ n2)2  X2(  X2  ♦ n2)3 


Ca 


6 


SlALi2  ♦,  n2)2  Pn  X(  A*  j n *Y 

sCa10 


i2  . _2*3 


Mat  A2I  xy  n2)2  + MatD  X*(  X*  ♦ n2) 
la6  ICa 


2/  i2  . 2v3 


Vl  * " °’ 


(13) 
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in  which  X • mrra/7  , K "A  - B . and  in  which  2K_  is  to  be  aubsti- 

n n n u 

tuted  for  Kq  at  n ■ 1;  Kn  * 0 when  n > n^.  When  equations  (11)  and  (12) 
are  added,  we  have 


D(  X2  + n2)4 

8 


D \ 2/  \2  2\2 
Nx  A (A  + n ) 

g 


2x3 


KXD  X ( A + n ) 
Ca8 


. Tn  \ ( \2  + n2)2  TDn  \(  \2  + p2)3 

*Ca10 


K* 

n 


Hat  A2(  l2  ♦ a2)2  4 Mate  A2(  >?  ± n2)3 
la6  Ka8 


<K’n-l  * KW  ' °- 


(U) 


in  which  K*n  * + Bq.  Again  at  n ■ 1,  2K?q  is  to  be  substituted  for 
K*0  and  K*n  * 0 when  n > n^. 

The  substitution  of  n from  0 to  in  equations  (13)  and  (14)  re- 
sults in  2a^  ♦ 1 simultaneous  algebraic  equations  for  2n^  + 1 unknowns. 
One  solution  to  these  equations  is  of  course,  the  trivial  one,  nanely. 


Kn  " K’n  " 0 (n  " °»1»2»3»  • • • t ni)» 


The  n on-trivial  solution  is  carried  out  in  the  following  section. 


-10- 


De termination  of  the  Buckling  Loads 


To  find  the  non-trivial  solutions  of  equations  (13)  and  (14)  in 

their  present  forms  is  a very  difficult  problem.  However,  from  tne 

12  3 

previous  investigations  * * it  was  found  that  the  minimum  buckling 

g 

loads  occur  at  X * 00  • Dividing  equations  (13)  and  (14)  by  X , and 
remembering  that  n^  may  also  be  very  large,  we  shall  drop  all  terms  con- 
taining X to  a power  greater  than  zero  in  the  denominator  but  keep 
terms  containing  the  ratio  n^/X  • Equations  (13)  and  (14)  thus  become 


(1  * x2  ' o'  ’ K”  * (K"-1  * K”+l)  ’ ° 


(15) 


2 H, 


(1  ♦ ?5  - -1  ♦ -V  ?)*'.♦—  * I’ntl>  ■ 0 


X2  C mA  X n 


IC 


(16) 


After  dividing  these  equations  through  by  Mat/lC , we  have 


K , + L K + K _ ■ 0, 
n+1  n n n-1  ' 


K»  ♦ L’K'  ♦ K*  . - 0 

n+x  n n n— 1 


(17) 

(18) 


in  which 


Mat 


va 


(19) 


L'_  ■ 


IC 

Mat 


_2  N.  T 

(1 , s.  . -1  , -X_  n ) 


X‘ 


C va 


(20) 


Following  the  examples  worked  out  in  Reference  3,  one  may  attempt  to 
solve  equations  (17)  and  (18)  as  finite-difference  equations.  This,  how- 
ever, was  not  found  possible.  The  reason  is  that  equations  (17)  and  (18) 
are  now  finite-difference  equations  with  variable  coefficients  and  the 
solutions  of  such  equations  are  difficult  mathematical  problems.  In- 
stead we  shall  solve  the  problem  in  the  following  manner. 
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Let  us  first  Investigate  the  magnitude  of  n^.  If  n^  is  small  com- 
pared to  X then  n^/ X — 0 as  \ — >-  oo  and  equations  (17)  and  (18) 

reduce  to  the  governing  equations  in  the  case  of  combined  compression 
and  bending.  In  order  that  the  torsional  load  may  heve  any  effect  on 
the  buckling  phenomena,  n^  must  be  also  e large  number  so  that  the  ratio 

n^/X  may  remain  finite.  For  any  n close  to  n^,  say,  n « n^-1,  n^-2, 

o.o,  n-,-q,  (q  < < n, ),  it  is  obvious  that  such  terms  as  l/X  > 2/X  t 

• • « , q/ A vanish  as  X— >~Oo  . Thus,  for  the  equations  in  which  n ■ n.  , 

n^-1,  n , -2 , . e . , the  following  relations  holds 


1 

A 


) 


- L (21) 

Therefore,  if  in  the  series  for  w,  we  take  the  summation  of  terms 
from  n^~q  to  n^,  the  system  of  equations  will  be  as  follows* 


and 


LK  + K , - 0 

"i  v1 

K ♦ IX  i ♦ K_  , - 0 
n^  n^-i  n^-2 


K . IZ 
^-q+1  nx-q 


L'K*  + K'  , - C 

"l  V1 


(22) 


K*  ♦ L«K«  . ♦ K*  - 0 
n^  n^-1  n^-2 


(23) 


K»  . ♦ L»K» 
ru-q-»l  n, -q 

- j.  n 
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where  L*  «* 


22- 

Mat. 


2 « 
nl  N1 


X X A T nl  x 

^ n 2 ~ 2 "\ 

r c Tra  u X 


To  obtain  the  non-trivial  solution  of  (22)  we  set  the  determinant 
of  the  coefficients  Kn  equal  to  zero,  namely, 


L 1 0 0 0.  ....  0 0 

1 L 1 0 0 ....  0 0 

0 1 L 1 0 ....  0 0 


0 ....  -L 


Equation  (24)  can  be  rewritten  as 


- 0 


(24) 


L - 


L - 


L - . 


0 


(25) 


1 


L 


L 


By  taking  more  and  more  terms  in  the  expression  for  v,  we  have 
successively  larger  and  larger  continued  fractions  in  (25)  and  in 
each  case  the  roots  L are  found. 

The  question  new  arises  as  to  which  of  the  roots  obtained  in 
solving  equation  (25)  is  the  one  which  should  be  used  in  the  final 
analysis.  This  is  easily  determinable  from  relation  (21),  which,  when 
rewritten,  becomes 


h.-E- 

Lat 


(1  + 


X2 


(26) 
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0bvlou3ly  the  lowest  value  of  the  critical  bending  moment  M occurs  fox 
the  largest  value  of  therefore,  in  solving  equation  (25),  the  largest 
root  obtained  in  each  case  is  the  one  selected- 

It  is  noted  (Fig„  2)  that  the  largest  root  of  equation  (26)  approach- 
es the  v*lue  of  2„ 00  as  larger  end  larger  continued  fractions  are  considered, 
Therefore  the  value  L *=  2e00  is  taken  a3  the  solution,.  Equation  (26) 
therefore  becomes,  after  dividing  through  by  TC/Letr 


2Mat 

IC 


Tra^C 


In  order  to  find  the  lowest  possible  value  of  the  term  2Mat/lC, 
it  is  necessery  to  minimize  the  right-hand  side  of  (27)  with  respect- 
to  the  variable  n^/  X c It  may  be  pointed  out  that  n X actually 

takes  on  only  discrete  values  on  account  of  the  integral  character  of 
n^  and  of  A («  mrra/Z  )„  However,  where  X is  very  large  (as  it  is  in 

this  ease),  it  is  possible  to  consider  n./ X as  a continuous  function 

in  the  minimizing  process „ 

Having  established  that  it  is  permissible  to  consider  n^/ X as 
a continuous  variable,  let  us  minimize  the  right-hand  side  of  (27)  by 
setting  _ 


a(n1/X) 


va2C 


Hence 


7T82C 


T 

— T 

2ira  C 


Consequently,  equation  (27)  becomes 


Dl 

X -0, 


2 Mat 


C £ TTa  C 


OS. 


-H- 


or 


T_ 

2vbl*C 


+ 


/ _ Hi  _ 2»S± 

V c ic 


(31) 


Exactly  the  sane  result  is  obtained  if  we  solve  for  the  non- trivial 
solution  for  equation  (23). 

In  the  case  where  the  cylinder  is  under  axial  compression  alone, 

M • T ■ 0,  and  equation  (31)  becomes 


or 


- n/  i - hyc 


o, 


N, 


(32) 


Equation  (32)  is  exactly  the  one  obtained  in  Reference  1. 

In  the  case  where  the  cylinder  is  under  torsional  loads  alone, 
M ■ • 0 and  equation  (31)  becomes 


or 


T 

2rra2C 


N « + C. 
*7  * 


(33) 


Equation  (33)  is  the  identical  equation  obtained  in  Reference  2. 
In  the  case  where  the  cylinder  is  under  bending  moments  alone, 
equation  (31)  becomes 


or 


- V 1 - 2Mat/lC  - 0, 

.c. 

l 


(34) 


Equation  (34)  is  exactly  the  one  obtained 
If  we  define  the  stress  ratios  R^, 
ing  formulas « 


in  Reference 
according  to 


3. 

the  follow- 
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"c 


Critical  compre -tnemf 

Buckling  stress  under  compression  alone 


2'  t Critical  bending  moment 

IC  Buckling  moment  under  bending  alone 


®T 


Critical  torsional  moment 
Buckling  moment  under  torsion  alone 


then  equation  (31)  may  ba  written  as 

Rt  - - \/l  - (Rg  ♦ Rq).  (35) 

The  interrelationship  between  compression,  bending,  and  torsional 
stress  ratios  given  by  equation  (35)  is  plotted  for  engineering  use  in 
Figs.  3 and  4*  Once  any  two  stress  ratios  are  specified,  the  buckling 
value  of  the  remaining  stress < ratio  can  be  determined  graphically  from 
these  curves. 
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FIGURE  I.  SANDWICH  CYLINDER 
CONFIGURATION  AND  LOADING 


Rb+  Rc 

FIGURE  3.  INTERACTION  CURVE:  TORSION  STRESS 
RATIO  vs.  SUM  OF  BENDING  AND  COMPRESSION 


STRESS  RATIOS 


FIGURE  4,  INTERACTION  CURVES: 


COMBINED  LOADING 


